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of  the  almrtna  compact  at  any  point,  and  to  produce  report  ready  graphs 
depicting  the  relationship  between  any  two  prescribed  variables. 

Analysis  of  other  errors  associated  with  the  use  of  the  dllatometer  shows 
that  these  are  negligible  compared  with  thermal  expansion  effects. 


The  rate  of  denslflcation  is  controlled  by  an  Interface  reaction  mechanism 
never  previously  observed  In  the  denslflcation  of  alumina.  Mass  transport 
Is  limited  by  the  movement  of  grain  boundary  dislocations  which  act  as  sites 
for  atoms  to  detach  from  grains.  The  actual  rate  limiting  process  is  the 
diffusion  of  solute  In  the  lattice  since  the  motion  of  solute  atoms  can 
result  In  a  large  number  of  atoms  being  freed  from  a  grain  boundary  dislo¬ 
cation.  Once  separated  from  the  dislocation  the  atoms  quickly  diffuse  away. 

Good  correlation  between  experimental  data  and  the  theory  Is  given  In  the 
relative  density  range  0.64-0.95,  although  the  theory  does  not  adequately 
describe  pore  geometry  above  0.9  relative  density.  The  rate  of  denslflcation 
In  the  temperature  range  1000-1150oC  Is  given  by: 
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where  D  “relative  density 

D0  *  Initial  relative  density 

R  *  particle  radius 

t  ■  time 

bb  ■  Burgers  vector  . 

bn  *  component  of  Burgers  reactor  normal 
to  the  grain  boundary 
p*  *  effective  contact  stress 

Ci  =  0.5 

Ds  *  diffuslvlty  of  solute 

n  *  effective  atomic  volume  of  solute 

G  *  shear  modulus 

e  *  ratio  of  solute  concentration  near  dislocation 

to  solute  concentration  In  bulk 
k  ■  Boltzmann's  constant 

T  ■  temperature 

C0  *  bulk  solute  concentration  measured 


The  measured  pressure  dependence  of  2.1  Is^also  In  good  agreement  with  the 
theoretical  prediction  of  2.  The  activation  energy  was  determined  to 
be  290kJ/mole. 
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SUMMARY  OF  DATA  AND  DATA  REDUCTION 


In  our  HIP  experiments*  we  continuously  measure  temperature* 
pressure*  and  change  In  can  diameter.  Changes  In  can  diameter  are  moni¬ 
tored  using  two  tungsten  probes  (see  Figure  1)  which  form  part  of  a  LYOT 
circuit.  We  calculate  a  continuous  record  of  density  as  a  function  of 
time  from  these  data.  As  a  result*  we  are  able  to  collect  much  more 
Information  from  a  single  run  than  Is  possible  In  ordinary  HIP  experi¬ 
ments  In  which  density  Is  recorded  only  at  the  beginning  and  end.  In 
this  section  all  the  experimental  data  Is  presented.  Furthermore*  we 
describe  the  procedure  for  converting  records  of  temperature*  pressure 
and  can  shrinkage  Into  a  continuous  record  of  density  as  a  function  of 
time  during  the  HIP  experiment.  Corrections  for  thermal  expansion  together  with 
assumptions  regarding  the  geometry  of  the  sample  and  can  during  denslfl- 
catlon  are  discussed.  A  Fortran  computer  program  which  carries  out 
these  calculations  is  given  In  the  appendix. 

The  experiments  were  planned  as  one  or  more  periods  of  time 
during  which  the  sample  was  held  at  roughly  constant  temperature  and 
pressure.  A  summary  of  the  conditions  for  all  six  experiments  Is  given 
In  Table  1.  The  experiment  numbers  In  the  table  will  be  used  throughout 
this  report  to  distinguish  the  various  runs.  In  experiments  21*  22a* 
and  22c*  temperature  and/or  pressure  were  changed  during  the  experiment 
so  that  additional  kinetics  data  could  be  obtained. 

Tables  2  through  7  give  values  of  experimental  data.  The  data 
were  Initially  recorded  on  a  strip-chart  recorder.  Values  of  pressure* 
temperature*  and  can  shrinkage  were  read  from  the  charts*  correlated 
with  clock  time*  and  keyed  Into  a  computer  for  processing.  Temperature 
was  measured  at  the  sample  tray*  but  a  simple  thermal  calculation  shows 
that  the  temperature  at  the  center  of  the  sample  will  be  largely  equili¬ 
brated  with  the  temperature  at  the  surface  of  the  can  within  one  minute. 

"Can  shrinkage"  is  the  change  In  the  outside  diameter  of  the  can  as 
measured  by  the  probe.  Note  that  the  can  shrinkage  Is  sometimes  nega¬ 
tive  since  thermal  expansion  Is  measurable  at  temperatures  too  low  for 
denslflcatlon.  The  tables  also  Include  relative  density  as  calculated 
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TABLE  1.  EXPERIMENTAL  CONDITIONS  FOR  HIP  EXPERIMENTS 


Experiment 

Nunber 

Pressure 

(MPa) 

Temperature 

(C) 

14 

101 

1150 

19 

102 

1000 

20 

103 

1050 

21 

34 

1000 

70 

1000 

22a 

34 

1050 

100 

1050 

100 

1150 

22c 

36 

1050 

100 

1050 

100 

1150 

TABLE  2.  DATA  FOR  EXPERIMENT  14 


Clock 

time 

Elapsed 
time  (min) 

Pressure 

(MPa) 

Temperature 

(C) 

Can  shrinkage 
(mm) 

Relative 

density 

13:49 

0.0 

6.9 

38 

0.000 

0.6531 

14:02 

13.0 

20.7 

51 

0.000 

0.6535 

14:22 

33.0 

41.4 

67 

0.038 

0.6583 

14:31 

42.0 

44.1 

64 

0.044 

0.6589 

14:47 

58.0 

73.1 

64 

0.051 

0.6596 

14:56 

67.0 

94.5 

65 

0.056 

0.6602 

14:57 

68.0 

99.3 

65 

0.056 

0.6602 

15:40 

111.0 

93.9 

115 

-0.013 

0.6538 

16:40 

171.0 

98.9 

672 

-0.176 

0.6540 

17:06 

197.0 

99.3 

875 

-0.184 

0.6610 

17:20 

211.0 

99.6 

1000 

0.260 

0.7215 

17:26 

217.0 

100.0 

1050 

0.657 

0.7805 

17:33 

224.0 

100.3 

1100 

1.086 

0.8540 

17:39 

230.0 

100.7 

1150 

1.467 

0.9307 

17:50 

241.0 

100.7 

1150 

1.586 

0.9559 

18:09 

260.0 

100.7 

1150 

1.594 

0.9576 
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TABLE  3.  DATA  FOR  EXPERIMENT  19 


Clock 

time 

Elapsed 
time  (min) 

Pressure 

(MPa) 

Temperature 

(C) 

Can  shrinkage 
(mm) 

Relative 

density 

12:44 

0.0 

0.0 

29 

0.000 

0.6536 

13:00 

16.0 

8.3 

36 

0.000 

0.6538 

13:05 

21.0 

16.2 

42 

0.000 

0.6540 

13:25 

41.0 

30.3 

58 

0.032 

0.6581 

13:41 

57.0 

49.0 

59 

0.040 

0.6591 

14:08 

84.0 

102.7 

59 

0.056 

0.6609 

15:01 

137.0 

97.9 

285 

0.032 

0.6653 

15:09 

145.0 

97.9 

362 

0.000 

0.6642 

15:18 

154.0 

98.6 

445 

-0.032 

0.6633 

15:27 

163.0 

98.6 

529 

-0.063 

0.6626 

15:52 

188.0 

99.3 

763 

“0.127 

0.6641 

16:13 

209.0 

100.3 

946 

-0.127 

0.6716 

16:19 

215.0 

100.5 

993 

0.000 

0.6894 

16:20 

216.0 

100.7 

1000 

0.032 

0.6938 

16:25 

221.0 

100.7 

1004 

0.222 

0.7190 

16:40 

236.0 

100.8 

1003 

0.492 

0.7571 

16:50 

246.0 

101.4 

1002 

0.603 

0.7738 

17:00 

256.0 

101.4 

1004 

0.691 

0.7875 

17:06 

262.0 

101.4 

1000 

0.730 

0.7936 

17:29 

285.0 

101.7 

1002 

0.865 

0.8159 

17:45 

301.0 

101.9 

1003 

0.953 

0.8308 

18:00 

316.0 

102.0 

1003 

1.000 

0.8391 

18:17 

333.0 

102.0 

1003 

1.048 

0.8476 

18:30 

346.0 

102.0 

1005 

1.087 

0.8549 

18:50 

366.0 

102.4 

1008 

1.135 

0.8639 

18:55 

371.0 

102.4 

1007 

1.143 

0.8653 

19:00 

376.0 

102.4 

1006 

1.151 

0.8667 

19:05 

381.0 

102.4 

1006 

1.175 

0.8712 
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TABLE  4.  DATA  FOR  EXPERIMENT  20 


Clock 

time 

Elapsed 
time  (min) 

Pressure 

(MPa) 

Temperature 

(C) 

Can  shrinkage 
(mm) 

Relative 

density 

12:45 

0.0 

0.0 

26 

0.000 

0.6485 

12:57 

12.0 

5.5 

27 

0.000 

0.6485 

13:01 

16.0 

15.2 

36 

0.000 

0.6487 

13:12 

27.0 

20.7 

48 

0.000 

0.6491 

13:19 

34.0 

24.1 

53 

0.000 

0.6492 

13:36 

51.0 

41.4 

57 

0.032 

0.6529 

13:59 

74.0 

80.7 

58 

0.048 

0.6548 

14:02 

77.0 

84.8 

57 

0.056 

0.6557 

14:10 

85.0 

100.0 

57 

0.056 

0.6557 

14:47 

122.0 

97.2 

215 

-0.016 

0.6523 

15:02 

137.0 

97.9 

327 

-0.048 

0.6522 

15:24 

159.0 

98.4 

496 

-0.063 

0.6561 

15:57 

192.0 

99.3 

766 

-0.127 

0.6588 

16:09 

204.0 

100.0 

855 

-0.127 

0.6623 

16:17 

212.0 

100.3 

914 

-0.095 

0.6685 

16:20 

215.0 

100.5 

937 

0.000 

0.6810 

16:25 

220.0 

100.7 

973 

0.119 

0.6974 

16:29 

224.0 

101.0 

1000 

0.286 

0.7205 

16:35 

230.0 

102.7 

1050 

0.571 

0.7635 

16:40 

235.0 

102.7 

1054 

0.826 

0.8032 

16:45 

240.0 

102.7 

1050 

0.984 

0.8293 

16:50 

245.0 

102.7 

1053 

1.111 

0.8518 

16:55 

250.0 

102.7 

1054 

1.214 

0.8707 

17:00 

255.0 

102.7 

1054 

1.302 

0.8872 

17:05 

260.0 

102.7 

1053 

1.365 

0.8995 

17:10 

265.0 

102.7 

1053 

1.405 

0.9073 

17:15 

270.0 

102.7 

1055 

1.460 

0.9187 

17:20 

275.0 

102.7 

1054 

1.516 

0.9302 

17:25 

280.0 

102.7 

1054 

1.548 

0.9368 

17:30 

285.0 

102.7 

1053 

1.611 

0.9504 

17:40 

295.0 

102.7 

1054 

1.627 

0.9539 

17:50 

305.0 

102.7 

1052 

1.683 

0.9660 

17:55 

310.0 

102.7 

1052 

1.691 

0.9678 

TABLE  5.  DATA  FOR  EXPERIMENT  21 


Clock  Elapsed 
time  time  (min) 

Pressure 

(MPa) 

Temperature 

(C) 

Can  shrinkage 

(mm) 

Relative 

density 

10:09 

0.0 

0.0 

23 

0.000 

0.6547 

10:35 

26.0 

20.7 

45 

-0.024 

0.6527 

10:47 

38.0 

27.6 

51 

0.008 

0.6565 

10:51 

42.0 

31.0 

52 

0.008 

0.6565 

10:57 

48.0 

36.0 

53 

0.024 

0.6584 

11:24 

75.0 

34.5 

168 

0.000 

0.6591 

12:01 

112.0 

34.5 

432 

-0.079 

0.6586 

12:33 

144.0 

34.6 

686 

-0.159 

0.6587 

12:42 

153.0 

34.8 

753 

-0.175 

0.6594 

13:00 

171.0 

34.5 

888 

-0.210 

0.6607 

13:13 

184.0 

34.5 

992 

-0.214 

0.6645 

13:15 

186.0 

34.5 

1000 

-0.191 

0.6677 

13:20 

191.0 

34.5 

1000 

-0.159 

0.6715 

13:23 

194.0 

34.5 

1000 

-0.127 

0.6753 

13:28 

199.0 

34.5 

1000 

-0.095 

0.6791 

13:35 

206.0 

34.5 

1000 

-0.063 

0.6830 

13:40 

211.0 

34.5 

1000 

-0.032 

0.6870 

13:47 

218.0 

34.5 

1000 

0.000 

0.6909 

13:55 

226.0 

34.5 

1000 

0.032 

0.6950 

14:06 

237.0 

34.5 

1000 

0.063 

0.6990 

14:18 

249.0 

34.5 

1000 

0.095 

0.7031 

14:33 

264.0 

34.5 

1000 

0.127 

0.7072 

14:50 

281.0 

34.5 

1000 

0.159 

0.7114 

15:15 

306.0 

34.5 

1000 

0.191 

0.7156 

15:30 

321.0 

34.5 

1000 

0.206 

0.7177 

15:45 

336.0 

55.8 

984 

0.222 

0.7191 

15:48 

339.0 

71.0 

998 

0.222 

0.7197 

15:51 

342.0 

71.0 

984 

0.238 

0.7212 

15:57 

348.0 

66.2 

999 

0.262 

0.7252 

16:18 

369.0 

73.1 

991 

0.341 

0.7357 

16:20 

371.0 

71.7 

999 

0.349 

0.7372 

16:28 

379.0 

70.3 

999 

0.381 

0.7417 

16:37 

388.0 

70.3 

1000 

0.413 

0.7463 

16:48 

399.0 

70.3 

1002 

0.445 

0.7509 

16:56 

407.0 

70.3 

1002 

0.476 

0.7555 

17:09 

420.0 

70.3 

1001 

0.508 

0.7601 

17:22 

433.0 

70.3 

1001 

0.540 

0.7649 

17:36 

447.0 

70.3 

1002 

0.571 

0.7697 

*■  - ' — 
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TABLE  6.  DATA  FOR  EXPERIMENT  22a 


Clock 

time 

Elapsed 
tine  (min) 

Pressure 

(MPa) 

Temperature 

(C) 

Can  shrinkage 

(am) 

Relative 

density 

00:30 

0.0 

0.0 

20 

0.000 

0.6362 

00:40 

10.0 

0.0 

20 

0.016 

0.6380 

00:50 

20.0 

0.0 

20 

-0.016 

0.6345 

00:55 

25.0 

0.0 

20 

0.016 

0.6380 

01:30 

60.0 

0.0 

22 

0.016 

0.6380 

01:40 

70.0 

8.3 

29 

0.003 

0.6368 

01:50 

80.0 

15.4 

42 

0.013 

0.6383 

02:00 

90.0 

22.9 

49 

0.016 

0.6388 

02:10 

100.0 

30.0 

51 

0.016 

0.6389 

02:20 

110.0 

33.1 

62 

0.022 

0.6399 

02:30 

120.0 

33.1 

103 

0.016 

0.6404 

02:40 

130.0 

33.1 

174 

-0.003 

0.6404 

02:50 

140.0 

33.1 

253 

-0.022 

0.6407 

03:00 

150.0 

33.1 

332 

-0.041 

0.6411 

03:10 

160.0 

33.1 

411 

-0.073 

0.6401 

03:20 

170.0 

33.1 

481 

-0.083 

0.6414 

03:30 

180.0 

33.1 

549 

-0.108 

0.6409 

03:40 

190.0 

33.1 

619 

-0.133 

0.6405 

03:50 

200.0 

33.1 

697 

-0.143 

0.6423 

04:00 

210.0 

33.1 

765 

-0.168 

0.6420 

04:10 

220.0 

33.1 

827 

-0.181 

0.6430 

04:20 

230.0 

33.1 

909 

-0.206 

0.6433 

04:21 

231.0 

33.1 

918 

-0.206 

0.6438 

04:24 

234.0 

33.1 

939 

-0.204 

0.6448 

04:25 

235.0 

33.1 

946 

-0.203 

0.6452 

04:27 

237.0 

33.1 

961 

-0.194 

0.6468 

04:30 

240.0 

33.1 

980 

-0.181 

0.6491 

04:33 

243.0 

33.1 

999 

-0.166 

0.6517 

04:35 

245.0 

33.3 

1013 

-0.156 

0.6534 

04:36 

246.0 

33.4 

1020 

-0.140 

0.6556 

04:39 

249.0 

33.8 

1046 

-0.092 

0.6623 

04:40 

250.0 

33.8 

1049 

-0.076 

0.6643 

04:42 

252.0 

33.8 

1054 

-0.020 

0.6713 

04:45 

255.0 

33.8 

1054 

0.063 

0.6816 

04:48 

258.0 

33.8 

1053 

0.102 

0.6864 

04:50 

260.0 

34.0 

1048 

0.127 

0.6893 

04:51 

261.0 

34.1 

1046 

0.138 

0.6906 

04:53 

263.0 

34.1 

1050 

0.159 

0.6936 

04:55 

265.0 

34.1 

1050 

0.181 

0.6964 

05:00 

270.0 

34.1 

1049 

0.235 

0.7033 

05:05 

275.0 

34.1 

1050 

0.292 

0.7110 

05:10 

280.0 

34.1 

1053 

0.333 

0.7167 

05:11 

281.0 

34.1 

1053 

0.338 

0.7174 

05:15 

286.0 

34.1 

1052 

0.359 

0.7201 

Vfc  **J*  4\-'  * 
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TABLE  6.  Continued 


Clock 

time 

Elapsed 
time  (min) 

Pressure 

(MPa) 

Temperature 

(C) 

Can  shrinkage 

(mm) 

Relative 

density 

05:20 

290.0 

34.1 

1050 

0.394 

0.7248 

05:21 

291.0 

34.1 

1050 

0.399 

0.7255 

05:25 

295.0 

34.1 

1047 

0.419 

0.7282 

05:30 

300.0 

34.1 

1043 

0.445 

0.7315 

05:31 

301.0 

34.1 

1042 

0.447 

0.7318 

05:35 

305.0 

34.1 

1053 

0.457 

0.7338 

05:40 

310.0 

34.1 

1053 

0.489 

0.7383 

05:50 

320.0 

34.1 

1052 

0.527 

0.7437 

06:00 

330.0 

34.1 

1052 

0.562 

0.7488 

06:10 

340.0 

34.1 

1052 

0.594 

0.7535 

06:20 

350.0 

34.2 

1051 

0.622 

0.7577 

06:30 

360.0 

34.2 

1050 

0.648 

0.7614 

06:40 

370.0 

34.3 

1049 

0.673 

0.7652 

06:50 

380.0 

34.3 

1048 

0.689 

0.7676 

07:00 

390.0 

34.4 

1046 

0.714 

0.7713 

07:10 

400.0 

34.5 

1045 

0.733 

0.7742 

07:20 

410.0 

34.5 

1047 

0.746 

0.7763 

07:30 

420.0 

34.5 

1049 

0.775 

0.7809 

07:40 

430.0 

34.5 

1051 

0.791 

0.7835 

07:50 

440.0 

34.5 

1052 

0.806 

0.7861 

08:00 

450.0 

34.5 

1051 

0.816 

0.7875 

08:10 

460.0 

34.5 

1051 

0.835 

0.7906 

08:20 

470.0 

34.5 

1051 

0.841 

0.7916 

08:30 

480.0 

34.5 

1052 

0.854 

0.7937 

08:40 

490.0 

34.5 

1052 

0.870 

0.7963 

08:50 

500.0 

34.5 

1051 

0.873 

0.7967 

09:00 

510.0 

34.5 

1050 

0.886 

0.7987 

09:10 

520.0 

34.5 

1049 

0.892 

0.7997 

09:20 

530.0 

34.5 

1048 

0.905 

0.8017 

09:30 

540.0 

34.5 

1047 

0.908 

0.8022 

09:40 

550.0 

34.5 

1052 

0.918 

0.6041 

09:50 

560.0 

37.2 

1047 

0.908 

0.8022 

09:52 

562.0 

40.0 

1043 

0.908 

0.8019 

09:56 

566.0 

45.5 

1050 

0.908 

0.8024 

09:58 

568.0 

46.3 

1045 

0.908 

0.8021 

10:00 

570.0 

51.7 

1047 

0.908 

0.8022 

10:02 

572.0 

55.2 

1049 

0.911 

0.8027 

10:05 

575.0 

60.3 

1048 

0.914 

0.8033 

10:06 

576.0 

62.1 

1048 

0.918 

0.8039 

10:09 

579.0 

68.9 

1047 

0.930 

0.8057 

10:10 

580.0 

72.1 

1047 

0.933 

0.8064 

10:13 

583.0 

81.4 

1047 

0.949 

0.8089 

10:14 

584.0 

77.2 

1050 

0.954 

0.8100 

10:15 

S8S.0 

77.8 

1049 

0.959 

0.8107 

10:20 

590.0 

80.9 

1045 

0.994 

0.8164 

■> 
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TABLE  6.  Continued 


Clock 

time 

Elapsed 
time  (min) 

Pressure 

(MPa) 

Temperature 

(C) 

Can  shrinkage 

(mm) 

Relative 

density 

10:23 

593.0 

82.7 

1042 

0.999 

0.8171 

10:25 

595.0 

87.3 

1045 

1.003 

0.8180 

10:26 

596.0 

89.6 

1047 

1.010 

0.8192 

10:29 

599.0 

96.5 

1047 

1.029 

0.8225 

10:30 

60C.0 

99.3 

1049 

1.035 

0.8237 

10:31 

601.0 

102.0 

1050 

1.042 

0.8249 

10:35 

605.0 

100.2 

1051 

1.070 

0.8299 

10:37 

607.0 

99.3 

1051 

1.083 

0.8321 

10:40 

610.0 

98.6 

1051 

1.102 

0.8354 

10:43 

613.0 

97.9 

1050 

1.119 

0.8384 

10:45 

61S.0 

97.9 

1050 

1.130 

0.8404 

10:50 

620.0 

97.9 

1051 

1.149 

0.8439 

10:53 

623.0 

97.9 

1051 

1.159 

0.8456 

10:55 

625.0 

97.9 

1051 

1.165 

0.8468 

11:00 

630.0 

97.9 

1051 

1.197 

0.8525 

11:10 

640.0 

99.8 

1051 

1.222 

0.8572 

11:20 

650.0 

101.0 

1050 

1.251 

0.8624 

11:30 

660.0 

100.7 

1051 

1.286 

0.8691 

11:40 

670.0 

100.7 

1052 

1.311 

0.8740 

11:50 

680.0 

100.7 

1052 

1.327 

0.8770 

12:00 

690.0 

100.7 

1051 

1.349 

0.8812 

12:10 

700.0 

100.7 

1052 

1.362 

0.8838 

12:20 

710.0 

100.7 

1052 

1.384 

0.8881 

12:30 

720.0 

100.7 

1051 

1.407 

0.8924 

12:40 

730.0 

100.7 

1050 

1.416 

0.8942 

12:50 

740.0 

100.7 

1051 

1.435 

0.8981 

13:00 

750.0 

100.7 

1051 

1.445 

0.9000 

13:10 

760.0 

100.7 

1049 

1.454 

0.9018 

13:20 

770.0 

100.7 

1052 

1.464 

0.9039 

13:30 

780.0 

100.7 

1051 

1.486 

0.9084 

13:40 

790.0 

100.7 

1051 

1.505 

0.9123 

13:50 

800.0 

100.7 

1061 

1.511 

0.9144 

14:00 

810.0 

100.7 

1074 

1.514 

0.9161 

14:10 

820.0 

100.7 

1094 

1.521 

0.9190 

14:20 

830.0 

100.7 

1104 

1.546 

0.9252 

14:30 

840.0 

100.7 

1134 

1.562 

0.9311 

14:40 

850.0 

100.7 

1153 

1.616 

0.9444 

14:50 

860.0 

100.7 

1150 

1.654 

0.9526 

15:00 

870.0 

100.7 

1151 

1.676 

0.9576 

15:10 

880.0 

100.6 

1153 

1.699 

0.9626 

15:20 

890.0 

100.3 

1150 

1.708 

0.9647 

15:30 

900.0 

100.3 

1151 

1.714 

0.9663 

15:40 

910.0 

100.3 

1152 

1.721 

0.9678 

15:44 

914.0 

100.3 

1152 

1.721 

0.9676 
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TABLE  7.  DATA  FOR  EXPERIMENT  22c 


Clock 

time 

Elapsed 
time  (min) 

Pressure 

(MPa) 

Temperature 

(C) 

Can  shrinkage 

(mm) 

Relative 

density 

01:18 

0.0 

34.5 

40 

0.000 

0.6272 

01:27 

9.0 

34.5 

68 

0.003 

0.6282 

01:36 

18.0 

34.5 

139 

-0.013 

01:48 

30.0 

34.5 

236 

-0.036 

0.6288 

01:57 

39.0 

34.5 

310 

-0.055 

0.6290 

01:58 

40.0 

34.5 

317 

-0.057 

0.6289 

02:07 

49.0 

34.5 

387 

-0.079 

0.6288 

02:16 

58.0 

34.5 

451 

-0.095 

0.6290 

02:25 

67.0 

34.5 

520 

-0.119 

0.6287 

02:34 

76.0 

34.8 

589 

-0.135 

0.6293 

02:43 

85.0 

35.2 

658 

-0.156 

0.6293 

02:52 

94.0 

35.2 

727 

-0.178 

0.6294 

03:01 

103.0 

35.2 

788 

-0.197 

0.6295 

03:09 

111.0 

35.2 

846 

-0.213 

0.6299 

03:19 

121.0 

35.2 

917 

-0.218 

0.6319 

03:22 

124.0 

35.2 

938 

-0.217 

0.6328 

03:25 

127.0 

35.2 

966 

-0.212 

0.6343 

03:28 

130.0 

35.2 

988 

-0.191 

0.6374 

03:37 

139.0 

35.2 

1050 

-0.022 

0.6583 

03:45 

147.0 

35.2 

1053 

0.168 

0.6801 

03:54 

156.0 

35.2 

1053 

0.296 

0.6955 

04:03 

165.0 

35.2 

1051 

0.382 

0.7061 

04:12 

174.0 

35.2 

1050 

0.453 

0.7151 

04:21 

183.0 

35.5 

1050 

0.507 

0.7220 

04:30 

192.0 

35.9 

1051 

0.555 

0.7283 

04:39 

201.0 

35.9 

1050 

0.591 

Hwk1! 

04:48 

210.0 

35.9 

1050 

0.618 

0.7368 

04:57 

219.0 

35.9 

1053 

0.659 

0.7424 

05:06 

228.0 

35.9 

1050 

0.687 

0.7461 

05:15 

237.0 

35.9 

1051 

0.704 

0.7484 

05:24 

246.0 

35.9 

1050 

0.735 

0.7528 

05:33 

255.0 

35.9 

1050 

0.758 

0.7560 

05:42 

264.0 

35.9 

1051 

0.782 

0.7594 

05:51 

273.0 

35.9 

1051 

0.800 

06:00 

282.0 

35.9 

1052 

0.810 

0.7635 

06:09 

291.0 

35.9 

1051 

0.826 

0.7656 

06:18 

300.0 

35.9 

1050 

0.845 

0.7683 

06:27 

309.0 

35.9 

1050 

0.861 

0.7707 

06:36 

318.0 

35.9 

1050 

0.875 

0.7727 

06:45 

327.0 

35.9 

1050 

0.889 

0.7748 

06:54 

336.0 

35.9 

1050 

0.903 

0.7768 

07:03 

345.0 

35.9 

1050 

0.914 

0.7785 

07:12 

354.0 

35.9 

1050 

0.927 

07:18 

360.0 

35.9 

1049 

0.933 

0.7613 

TABLE  7.  Continued 


Clock 

tine 

Elapsed 
tine  (nln) 

Pressure 

(MPa) 

Temperature 

(C) 

Can  shrinkage 
(mm) 

Relative 

density 

07:27 

369.0 

48.3 

1050 

—  .i 

___ 

07:36 

378.0 

62.1 

1046 

— 

— 

07:45 

387.0 

73.1 

1045 

— 

— 

07:54 

396.0 

89.6 

1046 

— 

— 

08:01 

403.0 

104.1 

1050 

1.102 

0.8071 

06:10 

412.0 

102.0 

1049 

1.172 

0.8182 

08:19 

421.0 

101.4 

1050 

1.226 

0.8269 

08:28 

430.0 

100.7 

1050 

1.274 

0.8348 

08:37 

439.0 

99.3 

1050 

1.312 

0.6412 

08:46 

448.0 

99.3 

1051 

1.347 

0.8473 

08:55 

457.0 

98.9 

1052 

1.373 

0.8517 

09:04 

466.0 

98.6 

1053 

1.398 

0.8562 

09:13 

475.0 

98.3 

1052 

1.433 

0.8621 

09:22 

484.0 

102.0 

1051 

1.457 

0.8663 

09:31 

493.0 

101.4 

1050 

1.485 

0.8712 

09:40 

502.0 

100.7 

1050 

1.506 

0.8751 

09:48 

510.0 

100.7 

1051 

1.525 

0.8786 

09:57 

519.0 

100.7 

1123 

1.535 

0.8855 

10:00 

522.0 

100.7 

1147 

1.556 

0.8910 

10:00 

522.5 

100.7 

1150 

1.561 

0.8921 

10:10 

532.0 

100.7 

1150 

1.660 

0.9110 

10:19 

541.0 

100.0 

1150 

1.722 

0.9232 

10:28 

550.0 

99.6 

1150 

1.750 

0.9287 

10:37 

559.0 

99.3 

1150 

1.765 

0.9318 

10:46 

568.0 

99.3 

1150 

1.773 

0.9333 

10:55 

577.0 

99.3 

1153 

1.778 

0.9346 

11:04 

586.0 

103.4 

1149 

1.784 

0.9356 

11:10 

592.0 

103.4 

1150 

1.789 

0.9367 

11:13 

595.0 

102.0 

1150 

1.789 

0.9367 

11:22 

604.0 

101.4 

1150 

1.789 

0.9367 

11:31 

613.0 

101.4 

1152 

1.789 

0.9368 
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by  the  methods  described  below. 

Before  experiment  22c*  a  change  was  made  In  the  gas  supply 
lines  for  the  autoclave.  This  design  change  apparently  resulted  In 
stronger  gas  currents  within  the  autoclave  during  pressurization.  The 
gas  currents  disturbed  the  probe  so  that  reliable  measurements  of  the 
can  shrinkage  could  not  be  obtained  during  pressurization  from  35.8  to 
104.1  MPa.  As  a  result*  listings  of  can  shrinkage  and  relative  density 
are  missing  for  a  small  portion  of  the  Table  7. 

Figures  2  through  7  reproduce  In  graphical  form  the  tempera¬ 
tures*  pressures*  and  densities  In  the  tables.  As  expected*  major 
changes  In  temperature  or  pressure  cause  clearly  visible  changes  In  the 
denslflcatlon  rate. 

Conversion  of  can  shrinkage  to  current  relative  density  also 
requires  the  measurements  shown  In  Table  8.  Figure  8  shows  the  posi¬ 
tions  at  which  the  various  dimensions  are  measured.  The  Initial  rela¬ 
tive  density  Is  obtained  from  the  Initial  weight*  length*  and  diameter 
of  the  sample.  Except  for  can  diameter  and  length*  the  Initial  dimen¬ 
sions  must  be  measured  before  assembly  of  the  can.  Most  of  the  final 
dimensions  require  cutting  of  the  can  and  sample.  While  six  measure¬ 
ments  are  necessary  to  describe  the  Initial  geometry*  only  four  are 
necessary  after  HIP  since  the  gap  between  sample  and  can  Is  gone. 

Reduction  of  data  starts  with  taking  data  for  temperature* 
pressure*  and  can  shrinkage  from  strip-chart  records  and  keying  It  Into 
a  computer.  Data  points  are  chosen  so  that  the  discrete  points  will 
provide  a  good  description  of  the  actual  experiment.  Since  the  time 
Interval  between  the  data  points  Is  generally  small*  linear  Interpola¬ 
tion  Is  used  between  readings.  Physical  dimensions  and  Initial  density, 
as  shown  In  Table  8,are  Inserted  at  the  beginning  of  the  data  file. 

The  data  Is  converted  to  densities  by  a  Fortran  program  that 
Is  listed  In  the  appendix.  Before  discussing  the  algorithm  used  In  the 
program*  however*  It  Is  appropriate  to  describe  the  assumptions  made  and 
the  estimated  error  due  to  the  assumptions. 

The  Integrated  thermal  expansions  of  alumina  and  stainless 
steel  are  given  by 


TIM  Mnl 


I 

Figure  2.  Oat*  for  experiment  14.  a)  Preaaure  (solid  line)  and 
temperature  (dotted  line),  b)  Density  as  calculated  from 
temperature  and  can  shrinkage. 


Figure  5.  Data  for  experiment  21.  a)  Pressure  (solid  lino)  and 
temperature  (dottod  lino),  b)  Donslty  as  calculatod  from 
temperature  and  can  shrinkage.  Donslty  curve  clearly  shows 
effect  of  changing  pressure  front  34  to  70  MPa  at  330  min. 
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Flgura  6.  Data  for  sxpsrlmsnt  22a.  a)  Pressure  (solid  11ns)  and 
temperature  (dottsd  11ns).  b)  Dsnslty  as  calculatsd  from 
temperature  and  can  shrlnkags.  Dsnslty  curvs  clsarly  shows 
offset  of  changing  prsssurs  from  34  to  100  MPa  at  580  min 
and  changing  tsmpsraturs  from  1050  to  115O0C  at  830  min. 


TIm  tain) 


Figure  7.  Data  for  experiment  22c.  a)  Pressure  (solid  lino)  and 
temperature  (dotted  lino),  b)  Density  as  calculated  from 
temperature  and  can  shrinkage.  Density  curve  clearly  shoes 
effect  of  changing  pressure  from  36  to  100  MPa  at  390  min 
and  changing  temperature  from  1050  to  1150°C  at  520  min. 
Density  Is  not  totally  reliable  between  361  and  396  minutes 
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TABLE  8.  PHYSICAL  MEASUREMENTS  OF  HIP  SAMPLES  AND  CANS 


Measurement*  mm 

(except  density) 

Run 

number 

14 

19 

20 

21 

22a 

22c 

Initial  sample  length 

39.67 

39.68 

39.70 

39.66 

39.67 

39.67 

Initial  can  length 

51.32 

52.24 

49.90 

50.91 

50.82 

50.95 

Final  can  length 

46.82 

47.90 

45.28 

49.04 

46.35 

46.50 

Initial  end  plug  length 

11.10 

11.13 

9.45 

10.97 

11.10 

11.10 

Final  end  plug  length 

11.69 

11.60 

10.08 

11.16 

11.79 

11.64 

Initial  sample  diameter 

22.00 

22.00 

22.00 

22.00 

22.01 

22.01 

Initial  can  diameter 

25.39 

25.37 

25.38 

25.36 

25.37 

25.37 

Final  sample  diameter 

23.28 

23.85 

23.33 

24.53 

23.11 

22.88 

Initial  can  wall  thickness 

1.57 

1.57 

1.57 

1.57 

1.66 

1.66 

Final  can  wall  thickness 

2.12 

1.98 

2.08 

1.83 

2.18 

2.11 

Initial  relative  density 

.6562 

.6556 

.6556 

.6549 

.6370 

.6374 

sample  (shaded  area)  is  exaggerated. 
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eAl2a5«  7.2978E-6  (T-293)  +  7.8486E-10  (T-293)  (T-800)  (la) 

e304ss  *  1.7751E-5  (T-293)  +  5.0423E-9  (T-293)  (T-800).  (lb) 

In  these  equations,  we  have  used  E  to  denote  multiplication  by  the 
specified  power  of  10.  following  the  standard  computational  convention. 
Note  that  these  are  are  the  Integrals  of  the  ordinary  coefficients  of 
linear  thermal  expansion.  These  equations  were  obtained  by  fitting  a 
quadratic  to  data  compiled  by  Toulouk1an(1) . 

Thermal  expansion  will  also  affect  the  geometry  of  the  probe 
mechanism  Itself.  Quite  a  few  thermal  effects  might  be  anticipated. 

The  probe  rods  will  change  In  length  and  diameter  both  above  and  below 
the  pivots.  The  pivots  themselves  will  expand  and  move.  The  sample 
will  be  displaced  as  Its  supporting  structure  expands.  Mhere  the  tem¬ 
perature  gradient  has  a  component  normal  to  the  axis  of  the  probe  rods, 
the  rods  will  deflect.  It  Is  not  feasible  to  perform  .  accurate  calcu¬ 
lation  of  all  the  thermal  expansion  effects,  so  It  would  appear  that  the 
best  approach  to  this  problem  would  be  to  calibrate  the  probe  carefully 
against  a  dummy  sample  of  some  non-denslflable  material  with  a  well- 
known  coefficient  of  thermal  expansion.  Even  if  such  a  calibration  were 
performed,  there  might  still  be  unconsidered  effects  from  thermal  tran¬ 
sients.  Our  calculations  to  date  have  used  the  simple  approach  of 
assuming  that  all  these  thermal  expansions  have  a  negligible  net  effect 
on  the  measured  can  diameter.  This  assumption  Is  supported  by  the 
apparently  small  systematic  errors  In  the  final  calculated  densities 
when  compared  to  Immersion  densities,  as  discussed  below. 

It  Is  assumed  that  the  material  of  the  can  Is  always  at  full 
density  and  that  the  mass  of  the  sidewall  of  the  can  Is  constant.  The 
stainless  steel  can  Is  rather  soft  at  HIP  temperatures,  so  Incompatibil¬ 
ity  of  thermal  expansion  between  sample  and  can  Is  handled  as  follows: 
For  any  given  density,  the  alumina  Is  assumed  to  expand  Isotropically 
upon  heating.  The  can  expansion  cannot  be  Isotropic,  however*  since  It 
expands  faster  than  the  alumina,  and  the  applied  gas  pressure  will  keep 
the  can  In  contact  with  the  sample.  Our  assumption  that  the  mass  of  the 
sidewall  Is  constant  Is  equivalent  to  assuming  that  the  can  does  not 
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slip  axially  along  the  sample,  and  that  the  extra  thermal  expansion  of 
the  can  Is  taken  up  as  a  plastic  thickening  of  the  can  wall.  Using  the 
final  dimensions  from  can  14*  the  combined  thermal  expansion  of  the  sam¬ 
ple  and  can  upon  heating  from  20°C  to  1150°C  will  be  0.378  mm.  This 
corresponds  to  a  change  In  the  relative  density  of  0.08*  which  Is  typi¬ 
cal.  Thermal  expansion  Is  undoubtedly  the  most  Important  single  correc¬ 
tion  In  the  conversion  of  probe  readings  to  densities. 

For  a  fixed  density  (expressed  as  fraction  of  theoretical  den¬ 
sity)  It  Is  assumed  that  the  thermal  expansion  of  alumina  Is  Isotropic. 
However*  experimental  measurements  show  that  the  length-to-dlameter 
ratio  of  the  sample  changes  as  the  density  changes.  Diameter  Is 
recorded  continuously*  but  length  Is  measured  only  at  the  beginning  and 
end  of  the  run.  Therefore*  It  Is  necessary  to  assume  some  relationship 
between  the  length  and  diameter  during  the  run.  It  Is  currently  assumed 
that 


where  d(p>  20)  and  L(p*  20)  denote  the  sample  diameter  and  length* 
respectively*  at  the  fractional  density  In  question  and  20°C.  The  sub¬ 
script  1  denotes  an  Initial  value.  The  exponent  a  Is  set  by  the  Initial 
and  final  measured  lengths  and  diameters  of  the  sample.  A  typical  value 
of  a  is  0.85.  This  Is  not  far  from  the  Ideal  value  a  *  1.0,  which 
describes  a  constant  length-to-dlameter  ratio. 

The  assumption  of  Equation  2  Is  essentially  arbitrary*  and 
other  relations  such  as 

L{p ,20)  -  L(pi ,20)  dfo,20)  -  d(pjt20) 

L(Pi,20)  -  L(Pf>  2d)  =  d(pf,20)  -  d(pr20)  (3) 

might  be  suggested*  where  the  subscript  f  denotes  a  final  value.  Using 
the  measurements  from  experiment  14*  one  finds  that*  for  any  diameter 
between  the  Initial  and  final  diameter*  the  maximum  discrepancy  In  cal¬ 
culated  sample  length  as  given  by  Equations  2  and  3  Is  about  0.013  mm. 
This  corresponds  to  an  error  In  density  of  less  than  0.041,  which  Is 
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negligible  In  our  experiments. 

It  Is  assumed  that  the  Initial  gap  between  the  sample  and  the 
can  shown  In  Figure  8  vanishes  during  HIP.  Thus  the  final  sample  diame¬ 
ter  Is  the  final  outside  diameter  of  the  can  minus  twice  the  final  can 
wall  thickness,  and  the  final  sample  length  Is  the  final  can  length 
minus  the  length  of  the  end  plugs.  This  assumption  Is  supported  by 
examination  of  the  samples;  no  gaps  were  found  after  HIP.  A  thin  reac¬ 
tion  layer  consisting  of  oxides  of  aluminum.  Iron,  chromium,  and  nickel 
was  detected  between  the  alumina  sample  and  the  stainless  steel  can. 
Measurements  of  the  reaction  layer  Indicate  that  an  error  of  less  than 
0.05%  results  from  Ignoring  the  layer. 

Using  these  assumptions,  the  program  proceeds  In  the  following 
way.  For  a  given  time,  the  pressure,  temperature,  and  shrinkage  are 
measured,  and  the  actual  outside  diameter  of  the  can  is  calculated  by 
subtracting  the  shrinkage  from  the  original  size.  Then  the  problem  Is 
attacked  from  the  other  direction:  A  fractional  density  Is  assumed,  and 
the  diameter  and  length  of  the  sample  at  20°C  are  calculated  by  means  of 
Equation  2.  Using  Equation  la,  the  dimensions  of  the  sample  at  tempera¬ 
ture  are  calculated.  The  volume  of  the  can  sidewall  at  20°C  Is 
corrected  for  temperature  by  converting  the  linear  expansion  of  Equation 
lb  to  a  volume  expansion,  and  the  necessary  wall  thickness  to  produce 
that  volume  Is  calculated.  For  the  assumed  density,  the  outside  diame¬ 
ter  of  the  can  will  be  the  sample  diameter  plus  twice  the  wall  thick¬ 
ness.  The  program  repeats  this  process  with  different  values  of  density 
until  It  finds  an  outside  diameter  which  Is  equal  to  the  actual  outside 
diameter  of  the  can.  As  output,  the  program  produces  a  table  In  which 
each  line  contains  time,  temperature,  pressure,  can  shrinkage,  and  rela¬ 
tive  density.  Plots  like  those  of  Figures  2  to  7  are  readily  made  by 
passing  selected  columns  to  a  plotting  package. 

One  difficulty  with  this  approach  Is  that  different  methods 
for  measuring  the  volume  of  the  can  sidewall  give  different  volumes. 

The  volume  of  the  sidewall  at  20°C  Is  given  by 


V  ■  irw(D-w)L 


(4) 


25 


where  w  Is  the  sidewall  thickness.  D  Is  the  outside  diameter  of  the  can. 
and  L  Is  the  length  of  the  sidewall,  all  Measured  at  20°C.  It  Is  rea- 
sonable  to  assume  that  the  volume  of  the  sidewall  as  measured  at  20°C  Is 
constant,  that  Is.  that  the  mass  of  the  sidewall  Is  constant.  The 
sidewall  volume  may  be  logically  defined  In  at  least  three  ways:  1.  w 
may  be  assumed  to  be  the  Initial  can  wall  thickness.  D  may  be  assumed  t# 
be  the  Initial  outside  diameter  of  the  can.  and  L  may  be  assumed  to  be 
the  Initial  can  length  minus  the  Initial  length  of  the  end  plugs.  2.  w 
may  be  assumed  to  be  the  Initial  can  wall  thickness.  D  may  be  assumed  to 
be  the  Initial  outside  diameter  of  the  can.  and  L  may  be  assumed  to  be 
the  Initial  sample  length.  3.  w  may  be  assumed  to  be  the  final  can  wall 
thickness.  D  may  be  assuaied  to  be  the  final  outside  diameter  of  the  can. 
and  L  may  be  assumed  to  be  the  final  sample  length.  Note  that  the 
length  used  In  method  1  Is  longer  than  the  length  used  In  method  2.  the 
difference  being  the  length  of  the  air  gaps  at  the  end  of  the  sample. 
Method  3  differs  from  the  others  by  using  final  measurements  Instead  of 
Initial  measurements.  Table  9  shows  the  sidewall  volumes  as  determined 
by  each  of  the  three  methods  Just  described.  As  expected,  the  volumes 
determined  by  method  2  are  consistently  smaller  than  those  determined  by 
method  1.  since  method  2  assumes  a  slightly  longer  sidewall.  However, 
there  Is  no  apparent  relationship  between  the  volumes  obtained  by  method 
3  and  those  obtained  by  methods  1  or  2.  The  discrepancies  between  the 
various  volumes  are  unexpectedly  large.  The  choice  of  the  sidewall 
volume  affects  the  density  that  Is  calculated  from  the  can  shrinkage. 

For  the  data  obtained  to  date,  the  three  methods  for  measuring  sidewall 
volume  can  yield  calculated  relative  densities  that  differ  by  as  much  as 
4.6X.  The  choice  of  the  correct  can  sidewall  volume  Is  an  Important 
unsolved  problem.  It  must  be  studied  In  more  detail  If  the  probe  Is  to 
provide  more  accurate  values  of  density. 

We  have  used  method  3  to  obtain  the  sidewall  volume,  since 
this  method  requires  no  assumptions  concerning  the  nature  of  the  col¬ 
lapse  of  the  can  onto  the  sample.  As  shown  In  Table  10.  there  has  gen¬ 
erally  been  good  correlation  between  the  final  density  as  calculated 
from  the  probe  data  and  as  measured  by  immersion.  If  experiment  22c  Is 
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TABLE  9.  CAN  SIDEWALL  VOLUMES  AS  DETERMINED 
BY  VARIOUS  METHODS 


Sidewall  volume,  cm^ 

Run 

Method  1 

Method  2 

Method  3 

14 

4.73920 

4.67457 

4.96023 

19 

4.83961 

4.67128 

4.94096 

20 

4.76382 

4.67497 

4.88847 

21 

4.69937 

4.66738 

4.95432 

22a 

4.92189 

4.91686 

4.95809 

22c 

4.93921 

4.91686 

4.80638 

I 


TABLE  10. 

COMPARISON  OF  FINAL  DENSITIES  AS 
CALCULATED  FROM  PROBE  DATA  AND 

AS  MEASURED  BY  IMMERSION 

Experiment 

Final  Density  Difference 

Number 

Calc. 

Maas. 

14 

.9576 

.960 

-.002 

19 

.8712 

.876 

-.005 

20 

.9678 

.952 

.016 

21 

.7697 

.758 

.012 

22a 

.9678 

.984 

-.016 

22c 

.9368  1 

.000 

-.063 
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excluded*  the  root-mea n- sq  ua  re  difference  between  the  calculated  and 
Immersion  densities  Is  1.3*.  There  Is  no  evident  systematic  error*  so 
our  assumptions  concerning  thermal  expansion  appear  to  be  justified. 

The  unusually  large  discrepancy  for  experiment  22c  has  not 
been  explained  completely.  It  Is  known  that  there  were  certain  Irregu¬ 
larities  In  the  experimental  procedure.  For  example*  a  reliable  cali¬ 
bration  of  the  probe  to  the  Initial  diameter  of  the  can  was  not  obtained 
due  to  equipment  difficulties  during  pressurization.  We  have  Included 
the  data  for  this  experiment*  however*  since  the  changes  In  density 
should  be  reliable  even  though  the  values  are  systematically  In  error. 
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THEORY  OF  DENSIFICATION 

In  this  section  a  theory  of  denslflcatlon  by  grain-boundary 
diffusion  with  Interface  reaction  control  Is  developed.  We  discuss  the 
role  of  grain-boundary  dislocations  In  Uniting  the  rate  of  grain¬ 
boundary  diffusion  and  review  the  geometry  of  Interparticle  necks  and 
the  effect  of  neck  geometry  on  the  effective  contact  stress.  We  present 
the  results  of  comparing  this  theory  of  denslflcatlon  with  experimental 
data  and  argue  that  Interface-reaction  controlled  grain-boundary  diffu¬ 
sion  Is  the  mechanism  that  best  explains  the  data.  Limitations  of  the 
theory  and  potential  for  future  development  are  discussed.  A  glossary* 
Table  11*  Is  provided  as  a  quick  reference  for  the  variables  used  In 
this  section. 

Denslflcatlon  by  grain-boundary  diffusion  has  been  thoroughly 
treated  In  the  literature.  The  effect  of  grain-boundary  dislocations  In 
limiting  grain-boundary  diffusion  has  also  been  discussed*  though  usu¬ 
ally  with  regard  to  creep  rather  than  denslflcatlon.  The  effect  of 
grain-boundary  dislocations  may  be  pictured  roughly  as  follows.  Con¬ 
sider  two  particles  In  a  sample  that  Is  denslfylng  by  grain-boundary 
diffusion.  Denslflcatlon  occurs  when  matter  flows  from  the  contact  area 
between  the  two  particles  to  the  surface  of  the  neck*  and  the  two  parti¬ 
cles  approach  each  other.  Under  certain  conditions*  the  rate  determin¬ 
ing  step  Is  not  diffusion  but  the  detaching  of  atoms  from  the  grains. 

The  energetically  favored  place  for  an  atom  to  detach  from  a  grain  Is  at 
a  grain-boundary  dislocation.  If  the  two  grains  were  separated.,  the 
grain-boundary  dislocations  would  be  ledges  on  the  new  surface.  As  den¬ 
slflcatlon  proceeds*  the  dislocation  climbs  along  the  grain  boundary. 

As  shown  In  Figure  9*  solute  atoms  tend  to  segregate  to  the  dislocation* 
and  climb  may  be  limited  by  solute  drag  If  the  dlffuslvlty  of  the  solute 
Is  low. 

The  terminology  "Interface-reaction  controlled"  is  actually  a 
misnomer.  While  It  is  true  that  the  amount  of  mass  transport  by  grain¬ 
boundary  diffusion  Is  limited  by  the  difficulty  of  moving  grain-boundary 
dislocations*  the  actual  rate-limiting  process  Is  the  diffusion  of 
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TABLE  11.  GLOSSARY  OF  TERMS 

average  neck  area 
Burgers  vector 
component  of  Burgers  vector  normal  to  grain  boundary 
derivative  of  coordination  number  with  respect  to  R»/R,  c  ■  15.5 
bulk  solute  concentration 
constant*  Cj  *  0.5 
relative  density 
Initial  relative  density 

relative  density  of  random  dense  packing*  DQ  *  0.64 

dlffuslvlty  of  solute  at  Infinite  temperature 

dlffuslvlty  of  solute 

activation  energy  for  diffusion  of  solute 

force  per  unit  length  on  dislocation 

shear  modulus 

Boltzmann's  constant 

dislocation  mobility 

external  pressure 

effective  contact  stress 

contact  pressure  due  to  surface  tension 

threshold  contact  pressure 
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TABLE  11.  Continued 

R  particle  radius 

R'  current  particle  radius  In  Arzt's  growing-sphere  description 
R  gas  constant 

T  temperature 

t  time 

vd1s  dislocation  velocity 

w  neck  axial  radius  of  curvature 

y  particle  center-to-center  distance 

Z  coordination  number 

ZQ  coordination  number  at  DQ,  ZQ  »  7.3 

o0  constant >  o0  *  0-2 

6  ratio  of  solute  concentration  near  dislocation  to  solute 

concentration  In  bulk 
y  surface  tension 

p  dislocation  density 

n  effective  atomic  volume  of  solute 


Figure  9.  Grain  boundary  dislocation  with  an  atmosphere  of  solute 
atoms.  As  the  dislocation  climbs  along  the  grain  boundary* 
solute  atoms  exert  a  drag  force. 
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solute  In  the  lattice.  The  motion  of  a  small  number  of  solute  atoms  In 
the  lattice  can  result  In  a  large  number  of  atoms  being  freed  from  a 
grain-boundary  dislocation.  Once  they  are  free  from  the  dislocation* 
the  atoms  quickly  diffuse  away. 

Following  Arzt*  Ashby*  and  Verrall^),  the  rate  at  which  an 
Individual  grain  boundary  dislocation  climbs  Is 

vd1s  «  MF  =  Mp«bn.  (5) 

Here  F  Is  the  force  per  unit  length  on  the  dislocation  and  M  Is  the 
dislocation  mobility.  From  elementary  dislocation  theory*  It  follows 
that  the  force  on  the  dislocation  Is  the  product  of  the  effective  con¬ 
tact  pressure  p*  and  the  component  of  the  Burgers  vector  normal  to  the 
grain  boundary  bR. 

Applying  the  results  of  Cottrell  and  Jaswon^3^  and  Cot- 
trell^*  Arrt*  Ashby*  and  Verrall  obtained  the  following  expression  for 
grain-boundary  dislocation  mobility  limited  by  solute  drags 

On  (6) 

y- 

ekTbfx 

b  o 

In  this  equation*  CQ  is  the  concentration  of  solute  In  the  lattice*  and 
the  concentration  In  the  solute  atmosphere  of  the  dislocation  Is 
Increased  by  a  factor  6  to  BC0.  Ds  Is  the  dlffuslvlty  of  the  solute  In 
the  lattice*  ft  Is  the  effective  atomic  volume  of  the  solute*  kT  has  Its 
usual  meaning*  and  b^  is  the  Burgers  vector  of  the  grain-boundary  dislo¬ 
cation. 

To  complete  the  description  of  a  grain-boundary  with  disloca¬ 
tions  we  need  only  an  expression  for  the  density  of  dislocations.  This 
problem  has  been  studied  by  Burton*5*  and  by  Arzt*  Ashby*  and  Ver- 
rall***.  Although  different  approaches  were  used*  their  results  are  In 
substantial  agreement.  Arzt's  result  Is 

cp.  m 

where  Cj  Is  a  constant  (about  0.5)*  6  Is  the  shear  modulus*  and  bb  Is 
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the  length  of  the  Burgers  vector. 

Before  proceeding  further*  It  Is  necessary  to  consider  again 
the  geometry  of  denslflcatlon  at  the  scale  of  complete  grains.  The  most 
Important  geometrical  quantities  are  Illustrated  In  Figure  10:  R  Is  the 
particle  radius*  y  Is  the  particle  center-to-center  distance*  x  is  the 
neck  radius*  and  tr  Is  the  neck  axial  radius  of  curvature.  If  It  is 
assumed  that  the  sample  starts  as  a  packing  of  tangential  spherical  par¬ 
ticles  and  that  the  particles  approach  each  other  as  denslflcatlon 
proceeds*  It  Is  clear  from  geometry  that 

D  -  D0(2R/y>3  (8) 

where  0o  is  the  Initial  density  and  0  Is  the  current  density.  By  dif¬ 
ferentiating  Equation  8*  one  obtains 


where  dy/dt  Is  the  rate  of  change  of  the  particle  center-to-center  dis¬ 
tance.  This  result  Is  at  variance  with  Equation  2  of  Arzt*  Ashby*  and 
Easterling*65.  Their  result  can  be  derived  from  the  assumption  that 

0  -  D0(y/2R)3.  (10) 

This  equation  Is  clearly  Incorrect;  It  predicts  a  decrease  In  density  as 
the  particles  move  toward  each  other.  Our  Equations  8  and  9  correct 
this  error. 

A  change  In  y  results  from  the  motion  of  grain-boundary  dislo¬ 
cations.  The  passage  of  one  dislocation  across  the  entire  neck  reduces 
the  center-to-center  distance  y  by  bn  where  bn  Is  the  component  of  the 
Burgers  vector  normal  to  the  grain  boundary.  In  the  case  where  several 
dislocations  are  moving  together*  the  number  of  dislocations  that  will 
pass  a  given  point  on  the  grain  boundary  In  unit  time  Is  v^sp.  Thus 
the  rate  of  approach  of  two  particles  Is 

dy/dt  ■  -  vd|s»bn  <u> 

It  has  been  Indicated  In  Equation  5  that  the  driving  force  for 


denslflcatlon  Is  the  effective  contact  pressure  p*  between  the  two  par¬ 
ticles.  Following  Arzt,  Ashby,  and  Easterl Ing*6*,  we  write  this  pres¬ 
sure  as 


n*  =  4”R  P  +  n  -  n 

P  aZD  ps  pth 


(12) 


where  p  Is  the  external  applied  pressure,  a  Is  the  average  neck  area,  Z 
Is  the  coordination  number,  ps  Is  the  contact  pressure  due  to  surface 
tension,  and  ptf>  Is  a  threshold  pressure  below  which  the  grain-boundary 
dislocations  cannot  move.  Since  our  material  was  thoroughly  outgassed 
before  canning  and  the  cans  welded  In  vacuum,  the  effect  of  gas  trapped 
In  the  pores  Is  neglected. 

Coordination  number  and  contact  area  have  been  treated  In  some 
detail  by  Arzt(7).  For  mathematical  convenience,  Arzt  describes  denslf¬ 
lcatlon  not  as  a  packing  of  spheres  that  approach  each  other  but  as  a 
packing  of  stationary  spheres  that  grow  fictitiously.  If  the  spheres 
grow  from  a  radius  R  to  a  radius  R',  the  Increase  In  the  volume  of  the 
spheres  will  cause  the  density  to  Increase  from  Its  Initial  value  DQ  to 
a  new  value  D,  where 


0/Do  «  (R*/R>3.  (l3) 

As  the  spheres  begin  to  grow,  It  Is  evident  that  they  must 
deform  to  keep  from  overlapping  with  their  neighbors.  Arzt  has  treated 
this  problem  In  the  following  way.  At  the  beginning  of  denslflcatlon, 
space  Is  partitioned  Into  Voronol  cells,  as  shown  In  Figure  11.  The 
Voronol  cell  of  a  given  particle  consists  of  all  points  that  are  closer 
to  the  center  of  that  particle  than  to  the  center  of  any  other  particle. 
Since  It  has  been  assumed  that  the  particles  are  stationary.  It  Is  also 
reasonable  to  assume  that  the  Voronol  cells  do  not  change  during  denslf¬ 
lcatlon.  Let  us  consider  a  single  sphere  as  shown  In  Figure  12.  The 
sphere  Is  tangent  to  the  boundary  of  its  Voronol  cell  at  any  point  where 
It  coordinates  another  sphere.  If  we  now  Imagine  the  sphere  growing 
without  deforming,  parts  of  the  sphere  will  fall  outside  the  Voronol 
cell,  as  shown  in  Figure  12b.  At  the  same  time,  we  must  also  Imagine 
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Figure  11.  A  two-dimensional  schematic  of  a  random  packing  of  spheres 
and  their  Voronol  cells*  after  Arzt. 


Figure  12.  A  typical  particle  In  a  random  dense  packing  and  Its 

Voronol  cell*  after  Arzt.  a)  No  denslflcatlon  has  occurred, 
b)  The  particle  grows  fictitiously  without  deformation  to 
radius  R' *  overlapping  neighboring  particles.  To  prevent 
overlap*  deformation  of  the  particle  occurs*  with  material 
being  distributed  either  c)  uniformly  over  the  free  surface 
or  d)  In  necks  around  the  contact  region. 
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the  coordinating  spheres  growing.  These  spheres  will  also  grow  beyond 
their  fixed  Voronol  cells*  and  the  spheres  will  overlap.  Real  particles 
do  not  overlap  during  denslflcatlon*  so  It  Is  clear  that  some  deforma¬ 
tion  of  the  spheres  must  take  place.  Arzt  has  assumed  that  when  a 
sphere  grows  to  the  point  that  part  of  Its  volume  would  be  outside  the 
Voronol  cell  boundary*  the  sphere  deforms  so  that  the  material  which 
would  fall  outside  the  cell  Is  redistributed  Inside  the  cell.  Under 
this  assumption*  one  may  think  of  the  sphere  as  growing  Inside  a  rigid 
Voronol  cell.  The  assumption  Is  actually  quite  reasonable.  Nherever 
sphere  A  would  grow  past  Its  cell  boundary  Into  the  cell  of  sphere  B* 
sphere  B  would  also  grow  past  Its  cell  boundary  Into  the  cell  of  sphere 
A.  Since  the  two  spheres  are  Identical*  they  would  force  each  other  to 
remain  within  their  respective  cells. 

Arzt  considered  two  ways  In  which  the  material  might  be  redis¬ 
tributed  when  the  spheres  deform.  Figure  12c  shows  the  material  uni¬ 
formly  spread  over  the  free  surface  of  the  particle.  This  would  be 
realistic  for  cold  compaction  of  a  ductile  powder.  Figure  12d  shows  the 
material  deposited  only  In  the  necks.  We  have  assumed  that  all  dis¬ 
placed  material  remains  In  the  necks*  since  It  Is  expected  that  surface 
tension  would  result  In  a  well-rounded  neck  like  that  shown  In  Figure 
12d  rather  than  a  sharp  neck  like  that  In  Figure  12c. 

As  the  spheres  continue  to  grow*  they  begin  to  coordinate 
additional  spheres  that  they  did  not  coordinate  In  the  original  packing. 
Using  radial  distribution  functions  for  a  random  dense  packing  developed 
by  Scott* 8)  and  Mason*®* »  Arzt  developed  the  following  approximate  equa¬ 
tion  for  the  coordination  number 

Z  -  Z0  +  c(RVR  -  1)  (14) 

where  Z0  -  7.3  and  c  »  15.5.  The  necks  which  result  from  the  original 
sphere  contact  points  will  be  larger  than  the  necks  which  form  later. 
Taking  this  fact  Into  account*  Arzt  has  derived  the  following  equation 
for  the  average  contact  area  If  all  material  remains  In  the  neckst 


a  -  11(Z0(R«  -  1)  +  c(R»  -  1)2/2)R2/ZR* 


(15) 
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With  this  Information*  the  first  term  of  Equation  12»  the  effect  of 
external  pressure*  may  be  evaluated* 

Following  Arzt*  Ashby*  and  Easterling*6* »  the  contact  pressure 
due  to  surface  tension  Is 

B  .  a  n  (16) 

ps  YU  X/ 

where  y  Is  the  surface  tension.  The  average  neck  area  Is  known  from 
Equation  15*  so  simple  geometry  will  give  the  neck  radius 

x  =  /ajt  (l7) 

A  reasonable  approximation  to  the  axial  radius  of  curvature  of  the  neck 
Is  given  by 

w  -  x2/(2R  -  2x) .  (18) 

This  equation  Is  not  strictly  consistent  with  Equation  IS.  Equation  15 
describes  a  neck  between  two  "overlapping”  spheres  like  the  one  shown  In 
Figure  13a:  the  material  In  the  neck  comes  from  the  volume  that  the 
spheres  would  overlap  If  they  were  complete.  Equation  18  describes  the 
simpler  geometry  of  Figure  13b  for  two  "tangent”  spheres:  the  material 
In  the  neck  Is  supplied  from  some  unspecified  source.  The  Inconsistency 
between  these  two  approaches  Is  Judged  to  be  Insignificant  for  the  fol¬ 
lowing  reasons:  First*  w  affects  only  p#,  which  will  be  shown  later  to 
be  of  secondary  importance  In  determining  p».  Second*  p#  Is  uncertain 
since  the  surface  tension  Is  not  well  known.  Third*  both  approaches 
assume  a  toroidal  neck.  Necks  of  this  type  ere  shown  In  the  axial 
cross-sections  of  Figure  13:  the  profile  of  the  neck  surface  Is  circular 
and  tangent  to  both  spheres.  A  toroidal  neck  Is  mathematically  con¬ 
venient*  but  this  geometry  would  not  be  expected  In  a  kinetic  system. 

For  these  reasons*  It  Is  concluded  that  the  expressions  for  w*  x»  and  p$ 
are  at  best  rather  rough  approximations. 

The  threshold  contact  pressure  for  motion  of  a  grain  boundary 
dislocation  Is  also  uncertain*  since  It  depends  on  the  Burgers  vector  of 
the  dislocations  and  the  type  of  dislocation  sources.  Ns  have  chosen  to 


use  the  approach  of  Arzt,  Ashby  and  Verrall.  The  threshold  stress  for 
grain  boundary  dislocation  notion  Is 


Gb,  (19) 

0  D 


where  aQ  =  0.2  Is  a  constant  which  reflects  the  relative  amount  that  the 
length  of  a  dislocation  must  change  In  order  to  move. 

As  Indicated  In  Equation  12,  the  effective  contact  pressure  Is 
the  sum  of  three  terms*  which  result  from  the  externally  applied  pres¬ 
sure*  surface  tension*  and  the  threshold  stress  for  grain-boundary 
dislocation  motion.  Let  us  use  the  following  data  to  calculate  the  mag¬ 
nitude  of  each  of  these  effects.  As  suggested  by  SEM  photographs  of  the 
fractured  HIP  specimens*  let  us  use  a  particle  radius  R  *  0.35E-6  m.  A 
surface  tension  y  *  1  J/m2  is  typical  for  ceramics.  The  data  of  Rysh- 
kewltch*10*  suggest  the  following  equation  for  the  shear  modulus: 

G  *  12.382E9  (1838.5-T)1/3  Pa  (20) 

for  temperatures  less  than  1600  K*  where  T  Is  the  temperature  In  kel- 
vlns.  We  have  used  a  Burgers  vector  for  lattice  dislocations  of  b  ■ 
4.75E-10  m  (for  slip  on  the  (0001)  <U20>  basal  system)*1**,  and,  fol¬ 
lowing  Arzt,  Ashby,  and  Verrall*2*,  let  bb  «  b/3.  Z0.  c,  and  <*0,  as 
specified  In  Table  II. 

Figure  14  shows  the  three  terms  of  the  effective  contact  pres¬ 
sure  as  functions  of  density.  These  curves  were  calculated  for  a  tem¬ 
perature  of  1150°C*  but  they  depend  only  weakly  on  temperature.  For  the 
sake  of  completeness,  the  curves  have  been  plotted  over  a  density  range 
from  0.65  to  0.995,  even  though  the  theory  Is  not  strictly  valid  over 
this  entire  range,  as  discussed  below.  The  effect  of  surface  tension  Is 
small  relative  to  the  effect  of  external  pressure  for  all  pressures  con¬ 
sidered.  At  low  densities,  the  magnitude  of  the  threshold  pressure  Is 
comparable  to  the  pressure  due  to  surface  tension.  At  a  density  of 
0.70,  an  external  applied  pressure  of  69  MPa  (10000  pel)  will  produce  a 
contact  pressure  of  507  MPa.  Surface  tension  effects  provide  an  addi¬ 
tional  30  IPs  of  contact  pressure,  and  the  threshold  pressure  Is  37  MPa, 
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Figure  14.  Contribution!  to  the  offoctlvo  contact  pressure  from 
thro*  sources  as  functions  of  density:  oxtsmal  prossuro 
(solid  linos) .  surface  tension  (dotted  line),  and  threshold 
for  grain-boundary  dislocation  motion  (dashed  line).  From 
bottom  to  top.  curves  for  external  pressure  are  for  34  MPa 
(5000  psl).  69  MPa  (10000  psl).  and  103  MPa  (15000  psl). 
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so  the  effective  contact  pressure  Is  500  MPa.  At  high  densities*  the 
magnitude  of  the  threshold  pressure  can  be  a  significant  fraction  of  the 
contact  pressure  due  to  the  applied  external  pressure.  Extending  the 
previous  example*  at  a  density  of  0.90*  an  external  applied  pressure  of 
69  MPa  will  produce  a  contact  pressure  of  99  MPa*  surface  tension 
effects  will  provide  only  2  MPa  of  contact  pressure*  and  the  threshold 
pressure  will  be  20  MPa*  producing  an  effective  contact  pressure  of  81 
MPa.  If  the  external  pressure  Is  cut  In  half*  the  effective  contact 
pressure  drops  to  32  MPa. 

The  total  effective  contact  pressure  due  to  all  three  terms  Is 
shown  In  Figure  15.  For  comparison*  the  bottom  curve  Is  drawn  with  no 
external  pressure*  corresponding  to  ordinary  sintering.  Note  that  the 
effective  contact  pressure  crosses  zero  at  a  density  near  0.69*  predict¬ 
ing  that  sintering  would  stop  at  this  density.  This  density  will  depend 
on  the  value  chosen  for  the  surface  tension  and  the  assumptions  concern¬ 
ing  neck  geometry  and  threshold  stress  for  dislocation  motion.  The 
theoretical  prediction  of  limited  sintering  Is  supported  by  our  work  on 
sample  preparation.  Before  HIP*  the  samples  were  sintered  for  3  hours 
at  1100°C.  The  relative  densities  obtained  were  typically  about  0.65; 
Table  8  lists  the  measured  Initial  densities. 

As  an  Initial  test*  It  Is  useful  to  ascertain  that  the  theory 
does  not  require  physically  unrealistic  values  for  any  of  Its  parame¬ 
ters.  Combining  Equations  5  to  7*  9*  and  11*  we  may  write 

do  aii/i)  yy-f  -•Vs0  <2I> 

"  IK\p0//  \bb )  Gb^SkTCo 

With  the  exception  of  Ds*  rough  estimates  can  be  made  of  all  of  the 
variables  In  this  equation.  Following  Arzt*  Ashby*  and  Verrall*  we  may 
take  bb  -  b/3»  bn  *  bb/VK  and  8  -  10.  The  value  of  n  is  unknown*  since 
It  Is  presently  not  known  what  solute  might  be  responsible  for  control¬ 
ling  the  motion  of  grain-boundary  dislocations.  However*  a  reasonable 
estimate  Is  n  ■  2.15E-29  m3,  which  is  half  of  the  volume  of  the  formula 
unit.  To  compare  theory  with  experiment*  let  us  duplicate  the  condi¬ 
tions  at  one  point  In  experiment  19.  At  an  elapsed  time  of  269  minutes* 
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we  find  that  0  ■  0.8,  T  *  1273  K,  p  ■  103  MPa,  and  the  denslftcatlon 
rate  dD/dt  *  1.6E-5  s”*.  When  all  these  values  are  substituted  Into 
Equation  21,  we  find  that,  for  the  theory  to  agree  with  the  experiment, 
we  must  have  Ds  *  1.2E-21  m*/s.  This  Is  not  an  unreasonable  dlf- 
fuslvlty,  so  for  this  one  piece  of  data.  It  has  been  shown  that  agree¬ 
ment  between  theory  and  experiment  can  be  achieved  with  physically  rea¬ 
sonable  values  for  the  parameters.  It  will  be  shown  below  that  these 
parameter  values  are  typical  of  those  obtained  by  curve-fitting. 

Although  the  controlling  solute  species  Is  not  known.  It  Is 
possible  to  measure  the  activation  energy  for  solute  diffusion,  E$.  By 
comparing  the  rates  of  denslflcatlon  at  1050°C  and  1150°C,  we  obtained 
the  values  Es  *  316  kJ/mol  for  experiment  22a  and  Es  «  263  kJ/mol  for 
experiment  22c,  or  an  average  value  of  290  kJ/mol .  Since  some  change  In 
density  occurred  while  the  temperature  was  rising,  a  better  value  for  Es 
would  be  obtained  If  the  measured  denslflcatlon  rates  were  corrected  for 
the  density  change.  Unfortunately,  a  meaningful  correction  could  not  be 
made  on  the  data  for  either  experiment.  In  the  case  of  experiment  22c, 
1150°C  was  not  reached  until  the  probe  Indicated  a  density  greater  than 
0.94.  This  density  Is  well  Into  the  final  stage,  and  we  have  not  yet 
extended  the  theoretical  model  to  the  final  stage.  In  experiment  22c, 
1150°C  was  reached  when  the  probe  Indicated  a  density  of  about  0.89, 
but,  as  discussed  in  the  previous  section,  the  calculated  densities  for 
this  experiment  are  probably  several  percent  below  their  actual  values. 
Since  the  actual  densities  are  not  known,  a  reasonable  correction  for 
effect  of  density  cannot  be  made.  Since  the  change  In  density  was 
small,  however,  the  density  correction  Is  expected  to  have  only  a  minor 
effect  on  the  calculated  activation  energy. 

An  additional  value  of  the  activation  energy,  Es  -  419  kJ/mol, 
was  obtained  by  comparing  the  data  for  experiments  19  and  20.  This 
value  Is  expected  to  be  less  accurate  than  those  given  above  since  It 
would  be  affected  by  any  variations  from  sample  to  sample. 

To  compare  the  predictions  of  theory  with  the  measurements  of 
experiment,  we  have  solved  the  differential  equation  for  denslflcatlon. 
Equation  21.  A  Runge-Kutta  method*1**  was  used  along  with  experimental 
records  of  external  pressure  and  temperature  as  listed  In  Tables  2  to  7. 
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Linear  Interpolation  was  used  between  tabulated  points.  For  the  pur¬ 
poses  of  Integration*  It  Is  convenient  to  recast  Equation  21  as  follows: 


dD  _  n/lV3  p*2exp(Es/RT) 
dt  ‘  \p0J  GT 


where  R  Is  the  gas  constant  and  n  Is  given  by 

,  /b. 


n  = 


n 


2R\b> 


2^l 

bbBkCo 


(23) 


In  Equations  22  and  23*  we  have  broken  down  the  solute  dlffuslvlty  using 
the  standard  Arrhenius  equation 

Ds  «  D°  exp(-Es/kT)  (24 

where  D°  Is  the  preexponent tal  factor  for  the  dlffuslvlty  of  the  solute 
and  R  Is  the  gas  constant.  Ne  have  defined  n  so  that  It  collects  the 
parameters  that  were  not  determined  by  our  experiments.  Except  for  n* 
all  the  variables  In  Equation  22  have  known  values:  the  Initial  density 
Dq,  the  temperature  T*  and  (through  Equation  12)  the  effective  contact 
pressure  p*  all  came  from  th'  -perlmental  data*  the  shear  modulus  G 
comes  from  the  literature*  and  the  activation  energy  for  solute  diffu¬ 
sion  Es  was  determined  by  experiment  as  discussed  above. 

We  have  used  n  as  the  single  adjustable  parameter  In  fitting 
Equation  22  to  our  data.  It  Is  appropriate  to  treat  n  as  adjustable* 
since  n  contains  several  variables  that  have  not  been  determined*  such 
as  bb  and  ft.  We  also  used  different  values  of  n  for  different  samples* 
so  the  theory  has  been  applied  to  each  experiment  Independently  rather 
than  to  all  six  experiments  as  a  group.  This  Is  also  appropriate  since 
n  contains  the  preexponentlal  factor  for  solute  diffusion  D|,  the  solute 
concentration  CQ*  and  the  ratio  of  dislocation  atmosphere  concentration 
to  bulk  concentration  0*  all  of  which  are  probably  subject  to  signifi¬ 
cant  sample-to-sample  variation. 

The  Initial  condition  for  Equation  22  Is 

D(t"0)  ■  D0. 
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This  was  the  Initial  condition  normally  used  In  our  Integrations.  In 
the  Integrations  for  experiments  22a  and  22c,  however.  Equation  25  could 
not  be  used.  When  the  density  D  *  Dr,  the  geometrical  description  of  * 

the  contacts  between  Individual  particles  breaks  down,  predicting  a  zero 
contact  area,  a  *  0.  Equation  12  then  gives  an  Infinite  effective  con¬ 
tact  pressure.  For  experiments  22a  and  22c,  DQ  <  Dr.  It  was  therefore 
necessary  to  choose  the  Initial  condition  so  that  D(t*0)  >  Dr,  even  If 
the  measured  Initial  density  was  lower  than  Df.  We  chose  to  Integrate 
with  an  Initial  density  of  0.641. 

The  choice  of  the  Initial  condition  Is  relatively  unimportant. 

This  can  be  seen  In  Figure  16.  Two  Integrations  were  performed  using 
the  temperature  and  pressure  data  for  experiment  14.  For  one  Integra¬ 
tion,  the  Initial  condition  was  as  specified  by  Equation  25.  For  the 
second  Integration,  we  arbitrarily  chose  D(tB0)  *  0.70.  The  same  value 
of  n  was  used  In  both  cases.  From  Figure  16  It  can  be  seen  that  even  a 
relatively  large  change  In  the  Initial  condition  rapidly  disappears  at 
denslflcatlon  proceeds.  The  reason  for  this  Is  straightforward.  If 
temperature  and  external  pressure  are  held  constant,  the  denslflcatlon 
rate  decreases  very  strongly  as  the  density  Increases  from  0.64  to  0.7. 

It  follows  that  a  change  In  the  Initial  density  will  largely  disappear 
when  significant  denslflcatlon  begins. 

Comparisons  between  data  and  experiment  are  shown  In  Figures 
17  to  22.  In  these  figures,  the  experimental  density  curves  of  Figures 
2b  to  7b  are  redisplayed  along  with  theoretical  density  curves.  The 
agreement  between  theory  and  experiment  Is  evident.  The  values  of  n 
used  for  all  of  these  calculations  are  relatively  close,  differing  by  at 
most  a  factor  of  6.  litis  Is  considered  to  be  acceptable  agreement. 

Since  n  Is  dependent  on  CQ,  any  variations  between  samples  In  the  con¬ 
centration  of  the  rate-controlling  solute  will  show  up  as  a  change  In  n. 

In  addition,  a  change  In  the  concentration  of  some  other  Impurity  may 
affect  D®  or  0,  and  thus  change  n. 

All  of  the  values  for  n  correspond  to  physically  reasonable 
choices  of  Dg,  bb,  ft,  etc.  This  can  be  seen  by  a  comparison  with  the 
sample  calculation  above.  As  has  been  argued,  the  values  of  the 
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Figure  16.  Effect  of  Initial  density  on  calculated  denelty  history. 
The  dotted  line  was  calculated  with  Initial  density  0.7, 
the  dashed  line  with  Initial  density  0.653.  For  both  curves 
n  ■  8E4.  The  temperature  and  pressure  histories  of 
experiment  14  were  used. 


■  r  •  |s4' 


Vjfc. 


A^TBUtQ 


Flflur*  17.  Comparison  between  experiment  and  theory  for  experiment 
14.  Solid  lino  Is  experiment'  dotted  lino  Is  theory. 
Theoretical  curve  was  calculated  with  n  •  8E4. 
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Figure  18.  Comparison  between  experiment  and  theory  for  experiment 
19.  Solid  line  Is  experiment*  dotted  line  Is  theory. 
Theoretical  curve  was  calculated  with  n  ■  3.2E4. 
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Figure  19.  Comparison  between  experiment  and  theory  for  experiment 
20.  Solid  line  Is  experiment*  dotted  line  Is  theory. 
Theoretical  curve  was  calculated  with  n  ■  8E4. 
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Figure  20.  Comparison  bwtween  experiment  and  theory  for  experiment 
21.  Solid  lino  Is  experiment*  dotted  lino  Is  theory. 
Theoretical  curve  was  calculated  with  n  •  1.4E4. 
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parameters  used  In  the  sample  calculation  are  reasonable*  and  they 
correspond  to  n  -  3.6E4.  To  provide  the  best  description  of  all  the 
data  of  experiment  19*  we  used  n  •  3.2E4.  In  either  case*  these  values 
of  n  are  typical  of  those  obtained  by  curve-fitting  for  all  of  the  six 
experiments*  so  the  values  of  n  used  In  the  Integrations  also  correspond 
to  physically  reasonable  choices  of  the  parameters. 

A  number  of  denslflcatlon  mechanisms  have  been  discussed  In 
the  literature*  Including  lattice  and  grain-boundary  diffusion  and  climb 
and  glide  of  dislocations.  None  of  these  mechanisms  describes  the 
results  of  our  experiments  as  well  as  does  the  theory  of  Interface- 
reaction  controlled  grain-boundary  diffusion  described  above.  A  summary 
of  the  characteristics  of  each  of  these  mechanisms  Is  shown  In  Table  12. 
Values  of  activation  energies  for  diffusion  come  from  Mohamed  and  Lang- 
don(13).  All  of  the  dlffuslonal  mechanisms  fall  In  that  they  do  not 
predict  the  correct  dependence  of  denslflcatlon  rate  on  pressure.  Stu¬ 
dies  of  data  for  experiments  21*  22a*  and  22c  Indicate  that  denslflca¬ 
tlon  rate  varies  like  p*  where  p  Is  the  external  pressure  and  the  pres¬ 
sure  exponent  m  Is  between  2  and  2.2.  By  contrast*  the  diffusion- 
controlled  mechanisms  all  predict  a  pressure  exponent  of  1.  The  dislo¬ 
cation  mechanisms  also  fall*  although  they  predict  a  pressure  dependence 
which  Is  stronger  than  that  of  Interface-controlled  diffusion.  There  Is 
also  no  notable  agreement  on  activation  energies*  although  the  energies 
for  grain-boundary  diffusion  are  moderately  close.  The  activation 
energy  for  dislocation  climb  Is  too  high.  The  activation  energy  for 
dislocation  glide  depends  on  the  slip  system*  but  at  the  temperatures 
used  In  our  experiments  (1000  to  1150°C)  only  the  basal  slip  system  Is 
active*  so  slip  could  contribute  only  a  small  amount  of  denslflcatlon. 
The  most  logical  explanation  of  the  data*  therefore*  Is  that  denslflca¬ 
tlon  occurs  by  Interface-reaction  controlled  grain-boundary  diffusion. 

As  It  presently  stands*  the  theory  described  above  would  be 
expected  to  break  down  outside  the  density  range  0.65  to  0.9.  At  densi¬ 
ties  below  Dr*  the  picture  of  a  random  dense  packing  is  Inadequate  to 
explain  the  density*  and  the  structure  of  an  aggregated  or  agglomerated 
powder  becomes  Important.  Denslflcatlon  of  an  agglomerated  powder  would 
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TABLE  12.  COMPARISON  OF  VARIOUS  MECHANISMS 
FOR  DENSIFICATION 


Mechanism 


Pressure  Activation  Energy 
Exponent  (kJ/mol) 


Lattice  Diffusion 

A1  Control 
0  Control 

Grain  Boundary  Diffusion 

A1  Control 
0  Control 

Dislocation  Climb 

Dislocation  Glide 

Interface-Reaction- 
Controlled  Diffusion 

Experiment 


1 

1 


1 

1 

-v  3 
5 

2 

2.1 


478 

637 

419 

226 

637 

? 

? 

290 


7 


be  difficult  to  model*  but  It  also  appears  that  It  Is  unnecessary  to 
model  It  since  HIP  would  tend  to  break  down  the  agglomerates.  When  the 
density  Is  only  slightly  above  Dr,  the  theory  Is  subject  to  doubt 
because  It  predicts  extremely  high  effective  contact  pressures.  At  high 
densities*  the  primary  difficulties  appear  to  be  geometrical.  Perhaps 
the  most  Important  effect  Is  that  the  pores  begin  to  close*  typically  at 
a  density  near  0.9*  and  a  completely  new  description  of  the  pore 
geometry  Is  needed.  It  should  be  noted*  however*  that  the  theory  Is  In 
good  agreement  with  experiment  to  densities  as  high  as  0.95.  From  an 
empirical  point  of  view*  the  model  Is  useful  In  the  density  range  0.9  to 
0.95  even  If  Its  theoretical  justification  Is  doubtful.  To  overcome 
these  difficulties*  we  currently  plan  to  extend  the  model  to  the  final 
stage  In  the  coming  year. 
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APPENDIX 


The  following  pages  contain  a  listing  of  a 
Fortran  program  to  convert  measurements  of 
can  shrinkage  Into  relative  densities, 
given  the  Initial  can  geometry.  The  program 
Includes  correction  for  thermal  expansion 
and  thickening  of  the  can  wall. 


o  o  o  o 
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program  probdens 
c 

c  This  program  converts  probe  readings  and  temperatures  to  sample 

c  densities, 

c 
c 

c  key  to  variable  names 

c  1  *  length 

c  d  =  diameter  (outside) 

c  w  *  wall  thickness 

c 

real  1  Inlt  sample.  1  flnl  sample.  1  Inlt  can, 

&  1  flnl  can,  1  Inlt  plugs,  1  flnl  plugs 

common  /sizes/  1  Inlt  sample,  1  flnl  sample,  1  Inlt  can, 

&  1  flnl  can,  1  Inlt  plugs,  1  flnl  plugs,  d  Inlt  sample, 

&  d  flnl  sample,  d  Inlt  can,  d  flnl  can,  w  Inlt  can, 

&  w  flnl  can,  density  Inlt,  alpha 

character*20  text 
c 

c  time  Is  In  minutes 

c  pressure  Is  In  psl 

c  temperature  Is  In  centigrade 

c  probe  reading  1$  In  (negative)  mils 

c  (positive  probe  reading  <«>  shrinkage) 

c  (negative  probe  reading  <«>  expansion) 

c 

c  dimensions  In  /sizes/  are  In  Inches  or  cubic  Inches 

c  alpha  and  density  Inlt  are  dimensionless 

c 

open  ( un1t**7 ,  f  11  e* '  canout * ,  status*  *  unknown  * ) 
open (un1t*5, status*' old') 
rew1nd(5) 

read  sample  dimensions 
call  readslze 

read  off  heading  lines 
read(5,'()«) 
read(5,'(a20)’)  text 
5  read(5, 50010, end*20)  time,  pressure,  temperature  c,  probe 

50010  format (5x,f7.1,f7.0,f6.0,f9.0) 

wr1te(6,6001)  time,  pressure,  temperature  c,  probe 
6001  format (5x,f5.1,f7.1,f6.1,f9. 2) 

temperature*temperature  c+273. 
c 

c  Correct  for  Initial  can  size. 

c  (d  can  Is  diameter  of  can  at  current  temperature.) 

c 


I 


d  can*d  Inlt  can-probe/1000. 

Find  upper  and  lower  limits  of  physical  probe  readings. 

density  min*. 6 
density  max*l. 

call  densprobCd  can  max, density  min, temperature) 
call  densprobCd  can  mlntdenslty  max, temperate re) 

1f(d  can.gt.d  can  max)  then 

wr1teC6,*C"  probe  reading  too  large.")') 
go  to  5 

else  If (d  can. It. d  can  min)  then 

wr1teC6,'C"  probe  reading  too  small.")') 
go  to  5 
end  If 


60010 


Calculate  density. 

1approx*0 

density  apx*dens1ty  m1n+(dens1ty  max-density  min)* 

( (d  can  max-d  can)/(d  can  max-d  can  min)) 
call  densprobCd  can  apx, density  apx, temperature) 

1apprax*1 approx+1 
If (d  can  apx.ge.d  can)  then 
d  can  max-d  can  apx 
density  m1n*dens1ty  apx 
else  IfCd  can  apx.le.d  can)  then 
d  can  m1n*d  can  apx 
density  max*dens1ty  apx 
end  If 

IfCCdensIty  max-density  mln.gt.l.e-5) .and. 1approx.lt. 12)  go  to  10 

wr1teC7, 60010)  time,  temperature  c,  pressure,  probe,  density  apx 

wr1teC6,60010)  time,  temperature  c,  pressure,  probe,  density  apx 

formatCf8.2,2f8.0,f8.2,f9.4) 

go  to  5 

stop 

end 

subroutine  densprobCd  can,dens1ty,temperature) 

This  routine  calculates  can  diameter  for  a  given  density. 

real  1  Inlt  sample,  1  flnl  sample,  1  Inlt  can, 

1  flnl  can,  1  Inlt  plugs,  1  flnl  plugs 

common  /sizes/  1  Inlt  sample,  1  flnl  sample,  1  Inlt  can, 

1  flnl  can,  1  Inlt  plugs,  1  flnl  plugs,  d  Inlt  sample, 
d  flnl  sample,  d  Inlt  can,  d  flnl  can,  v  Inlt  can, 
w  flnl  can,  density  Inlt,  alpha 
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real  1  sample#!  sample  Inlt#  1  sample  Inlt  20 
data  pl/3. 1415926/ 
c 

1  sample  Inlt  20-1  Inlt  sample 

volume  Inlt  20*p1#l  Inlt  sample»d  Inlt  sample*d  Inlt  sample/4. 

d  sample  Inlt  20-d  Inlt  sample 

vol  can  20-(l  flnl  can  -  1  flnl  plugs)*p1* 

4  (d  flnl  sample  +  w  flnl  can)*w  flnl  can 

c 

c  Calculate  Integrated  thermal  expansion, 

c 

expand  a12o3-7.2978e-6*(temperature-293.)+ 

&  7. 8486e-10« (temperature-293 . )*(temperature-800. ) 

expand  304ss-l .7751e-5# (temperature-293 . )+ 

&  5 .0423e-9»( temperature-293 . )* (temperature-600. ) 

c 

c  Correct  lengths  and  volumes  for  temperature, 

c 

d  sample  Inlt-d  sample  Inlt  20*(l.+expand  al2o3) 

1  sample  lnlt-1  sample  Inlt  20*(l.+expand  al2o3) 
volume  Inlt-volume  Inlt  20«(l.+3.#expand  al2o3) 
vol  can-vol  can  20*(l.+3.#expand  304ss) 
c 

c  Calculate  can  geometry, 

c 

volume-volume  1n1t*(dens1ty  Inlt/denslty) 

d  sample-((4«volume*d  sample  1n1t**alpha)/(p1»l  sample  1n1t))»» 
&  (l./(2.+alpha)> 

1  sample-(4.*volume)/(p1*d  sample**2> 

«a!1-(-p1*l  sampled  sample^ 

&  sqrt((p1*l  sample*d  sample>»»2*-4.,lp1»l  sample*vol  can)) 

&  /(2.*p1*l  sample) 

d  can*d  sample+2.*wal 1 
return 
end 

subroutine  readslze 
c 

c  this  routine  reads  size  dat  for  the  sample 

c 

real  1  Inlt  sample#  1  flnl  sample#  1  Inlt  can# 

A  1  flnl  can#  1  Inlt  plugs#  1  flnl  plugs 

common  /sizes/  1  Inlt  sample#  1  flnl  sample#  1  Inlt  can# 
t  1  flnl  can#  1  Inlt  plugs#  1  flnl  plugs#  d  Inlt  sample# 

A  d  flnl  sample#  d  Inlt  can#  d  flnl  can#  «  Inlt  can# 

4  v  flnl  can#  density  Inlt#  alpha 
c 

read  (5#505)  1  Inlt  sample#  1  Inlt  can# 

A  1  flnl  can#  1  Inlt  plugs#  1  flnl  plugs#  d  Inlt  sample# 


